Metagenomics for Pathogen Detection During a Mass Mortality Event in Songbirds.

Mass mortality events in wildlife can be indications of an emerging infectious disease. During the spring and summer of 2021, hundreds of dead passerines were reported across the eastern US. Birds exhibited a range of clinical signs including swollen conjunctiva, ocular discharge, ataxia, and nystagmus. As part of the diagnostic investigation, high-throughput metagenomic next-generation sequencing was performed across three molecular laboratories on samples from affected birds. Many potentially pathogenic microbes were detected, with bacteria forming the largest proportion; however, no singular agent was consistently identified, with many of the detected microbes also found in unaffected (control) birds and thus considered to be subclinical infections. Congruent results across laboratories have helped drive further investigation into alternative causes, including environmental contaminants and nutritional deficiencies. This work highlights the utility of metagenomic approaches in investigations of emerging diseases and provides a framework for future wildlife mortality events.

Endogenous feline leukemia virus long terminal repeat integration site diversity is highly variable in related and unrelated domestic cats.

Endogenous retroviruses (ERV) are indicators of vertebrate evolutionary history and play important roles as homeostatic regulators. ERV long terminal repeat (LTR) elements may act as cis-activating promoters or trans-activating enhancer elements modifying gene transcription distant from LTR insertion sites. We previously documented that endogenous feline leukemia virus (FeLV)-LTR copy number variation in individual cats tracks inversely with susceptibility to virulent FeLV disease. To evaluate FeLV-LTR insertion characteristics, we assessed enFeLV-LTR integration site diversity in 20 cats from three genetically distinct populations using a baited linker-mediated PCR approach. We documented 765 individual integration sites unequally represented among individuals. Only three LTR integration sites were shared among all individuals, while 412 sites were unique to a single individual. When primary fibroblast cultures were challenged with exogenous FeLV, we found significantly increased expression of both exogenous and endogenous FeLV orthologs, supporting previous findings of potential exFeLV-enFeLV interactions; however, viral challenge did not elicit transcriptional changes in genes associated with the vast majority of integration sites. This study assesses FeLV-LTR integration sites in individual animals, providing unique transposome genotypes. Further, we document substantial individual variation in LTR integration site locations, even in a highly inbred population, and provide a framework for understanding potential endogenous retroviral element position influence on host gene transcription.

Feline Leukemia Virus Frequently Spills Over from Domestic Cats to North American Pumas.

Feline leukemia virus (FeLV) is a gammaretrovirus with horizontally transmitted and endogenous forms. Domestic cats are the primary reservoir species, but FeLV outbreaks in endangered Florida panthers and Iberian lynxes have resulted in mortalities. To assess prevalence and interspecific/intraspecific transmission, we conducted an extensive survey and phylogenetic analysis of FeLV infection in free-ranging pumas (n = 641) and bobcats (n = 212) and shelter domestic cats (n = 304). Samples were collected from coincident habitats across the United States between 1985 and 2018. FeLV infection was detected in 3.12% of the puma samples, 0.47% of the bobcat samples, and 6.25% of the domestic cat samples analyzed. Puma prevalence varied by location, with Florida having the highest rate of infection. FeLV env sequences revealed variation among isolates, and we identified two distinct clades. Both progressive and regressive infections were identified in cats and pumas. Based on the time and location of sampling and phylogenetic analysis, we inferred 3 spillover events between domestic cats and pumas; 3 puma-to-puma transmissions in Florida were inferred. An additional 14 infections in pumas likely represented spillover events following contact with reservoir host domestic cat populations. Our data provide evidence that FeLV transmission from domestic cats to pumas occurs widely across the United States, and puma-to-puma transmission may occur in genetically and geographically constrained populations. IMPORTANCE Feline leukemia virus (FeLV) is a retrovirus that primarily affects domestic cats. Close interactions with domestic cats, including predation, can lead to the interspecific transmission of the virus to pumas, bobcats, or other feline species. Some infected individuals develop progressive infections, which are associated with clinical signs of disease and can result in mortality. Therefore, outbreaks of FeLV in wildlife, including the North American puma and the endangered Florida panther, are of high conservation concern. This work provides a greater understanding of the dynamics of the transmission of FeLV between domestic cats and wild felids and presents evidence of multiple spillover events and infections in all sampled populations. These findings highlight the concern for pathogen spillover from domestic animals to wildlife but also identify an opportunity to understand viral evolution following cross-species transmissions more broadly.

Rapid evolution of SARS-CoV-2 in domestic cats.

SARS-CoV-2 (SARS2) infection of a novel permissive host species can result in rapid viral evolution. Data suggest that felids are highly susceptible to SARS2 infection, and species-specific adaptation following human-to-felid transmission may occur. We employed experimental infection and analysis of publicly available SARS2 sequences to observe variant emergence and selection in domestic cats. Three cohorts of cats (N = 23) were inoculated with SARS-CoV-2 USA-WA1/2020 or infected via cat-to-cat contact transmission. Full viral genomes were recovered from RNA obtained from nasal washes 1-3 days post-infection and analyzed for within-host viral variants. We detected 118 unique variants at ≥3 per cent allele frequency in two technical replicates. Seventy of these (59 per cent) were nonsynonymous single nucleotide variants (SNVs); the remainder were synonymous SNVs or structural variants. On average, we observed twelve variants per cat, nearly 10-fold higher than what is commonly reported in human patients. We observed signatures of positive selection in the spike protein and the emergence of eleven within-host variants located at the same genomic positions as mutations in SARS2 variant lineages that have emerged during the pandemic. Fewer variants were noted in cats infected from contact with other cats and in cats exposed to lower doses of cultured inoculum. An analysis of ninety-three publicly available SARS2 consensus genomes recovered from naturally infected domestic cats reflected variant lineages circulating in the local human population at the time of sampling, illustrating that cats are susceptible to SARS2 variants that have emerged in humans, and suggesting human-to-felid transmission occurring in domestic settings is typically unidirectional. These experimental results underscore the rapidity of SARS2 adaptation in felid hosts, representing a theoretical potential origin for variant lineages in human populations. Further, cats should be considered susceptible hosts capable of shedding virus during infections occurring within households.

Prevalence and Genomic Sequence Analysis of Domestic Cat Hepadnavirus in the United States.

Hepadnaviruses are partially double-stranded DNA viruses that infect a variety of species. The prototypical virus in this family is the human hepatitis B virus, which chronically infects approximately 400 million people worldwide and is a risk factor for progressive liver disease and liver cancer. The first hepadnavirus isolated from carnivores was a domestic cat hepadnavirus (DCH), initially identified in Australia and subsequently detected in cats in Europe and Asia. As with all characterized hepadnaviruses so far, DCH infection has been associated with hepatic disease in its host. Prevalence of this infection in the United States has not been explored broadly. Thus, we utilized conventional and quantitative PCR to screen several populations of domestic cats to estimate DCH prevalence in the United States. We detected DCH DNA in 1 out of 496 animals (0.2%) in the U.S. cohort. In contrast, we detected circulating DCH DNA in 7 positive animals from a cohort of 67 domestic cats from Australia (10.4%), consistent with previous studies. The complete consensus genome of the U.S. DCH isolate was sequenced by Sanger sequencing with overlapping PCR products. An in-frame deletion of 157 bp was identified in the N-terminus of the core open reading frame. The deletion begins at the direct repeat 1 sequence (i.e., the 5' end of the expected double-stranded linear DNA form), consistent with covalently closed circular DNA resultant from illegitimate recombination described in other hepadnaviruses. Comparative genome sequence analysis indicated that the closest described relatives of the U.S. DCH isolate are those previously isolated in Italy. Motif analysis supports DCH using NTCP as an entry receptor, similar to human HBV. Our work indicates that chronic DCH prevalence in the U.S. is likely low compared to other countries.

Apathogenic proxies for transmission dynamics of a fatal virus.

Identifying drivers of transmission-especially of emerging pathogens-is a formidable challenge for proactive disease management efforts. While close social interactions can be associated with microbial sharing between individuals, and thereby imply dynamics important for transmission, such associations can be obscured by the influences of factors such as shared diets or environments. Directly-transmitted viral agents, specifically those that are rapidly evolving such as many RNA viruses, can allow for high-resolution inference of transmission, and therefore hold promise for elucidating not only which individuals transmit to each other, but also drivers of those transmission events. Here, we tested a novel approach in the Florida panther, which is affected by several directly-transmitted feline retroviruses. We first inferred the transmission network for an apathogenic, directly-transmitted retrovirus, feline immunodeficiency virus (FIV), and then used exponential random graph models to determine drivers structuring this network. We then evaluated the utility of these drivers in predicting transmission of the analogously transmitted, pathogenic agent, feline leukemia virus (FeLV), and compared FIV-based predictions of outbreak dynamics against empirical FeLV outbreak data. FIV transmission was primarily driven by panther age class and distances between panther home range centroids. FIV-based modeling predicted FeLV dynamics similarly to common modeling approaches, but with evidence that FIV-based predictions captured the spatial structuring of the observed FeLV outbreak. While FIV-based predictions of FeLV transmission performed only marginally better than standard approaches, our results highlight the value of proactively identifying drivers of transmission-even based on analogously-transmitted, apathogenic agents-in order to predict transmission of emerging infectious agents. The identification of underlying drivers of transmission, such as through our workflow here, therefore holds promise for improving predictions of pathogen transmission in novel host populations, and could provide new strategies for proactive pathogen management in human and animal systems.

Multiple Introductions of SARS-CoV-2 Alpha and Delta Variants into White-Tailed Deer in Pennsylvania.

The SARS-CoV-2 pandemic began by viral spillover from animals to humans; today multiple animal species are known to be susceptible to infection. White-tailed deer, Odocoileus virginianus, are infected in North America at substantial levels, and genomic data suggests that a variant in deer may have spilled back to humans. Here, we characterize SARS-CoV-2 in deer from Pennsylvania (PA) sampled during fall and winter 2021. Of 123 nasal swab samples analyzed by RT-qPCR, 20 (16.3%) were positive for SARS-CoV-2. Seven whole genome sequences were obtained, together with six more partial spike gene sequences. These annotated as alpha and delta variants, the first reported observations of these lineages in deer, documenting multiple new jumps from humans to deer. The alpha lineage persisted in deer after its displacement by delta in humans, and deer-derived alpha variants diverged significantly from those in humans, consistent with a distinctive evolutionary trajectory in deer. IMPORTANCE Coronaviruses have been documented to replicate in numerous species of vertebrates, and multiple spillovers of coronaviruses from animals into humans have founded human epidemics. The COVID-19 epidemic likely derived from a spillover of SARS-CoV-2 from bats into humans, possibly via an intermediate host. There are now several examples of SARS-CoV-2 jumping from humans into other mammals, including mink and deer, creating the potential for new animal reservoirs from which spillback into humans could occur. For this reason, data on formation of new animal reservoirs is of great importance for understanding possible sources of future infection. Here, we identify extensive infection in white-tailed deer in Pennsylvania, including what appear to be multiple independent transmissions. Data further suggests possible transmission among deer. These data thus help identify a potential new animal reservoir and provide background information relevant to its management.

Modelling the role of predation on disease burdens of prey.

Research Highlight: Brandell, E. E., Cross, P. C., Smith, D. W., Rogers, W., Galloway, N. L., MacNulty, D. R., Stahler, D. R., Treanor, J. & Hudson, P. J. (2022). Examination of the interaction between age-specific predation and chronic disease in the Greater Yellowstone Ecosystem. Journal of Animal Ecology, 00, 1-12. https://doi.org/10.1111/1365-2656.13661. Predation can alter disease dynamics in prey. If predators select for infected individuals, they can reduce disease burdens. In other cases, predators can increase disease burdens via various mechanisms such as altered prey behaviour. The influence of predation on disease dynamics is a result of interactions among various traits of the predators, prey and the pathogen itself. For example, pathogens tend to vary with age and predators typically select for certain age classes. Thus, the overlap between ages selected by predators and those infected will likely contribute to any effects of predation on reducing disease burdens. In this paper, Brandell et al. (2022) develop a model to evaluate the predator cleansing effect given age-based variation in pathogens and predation. The model was developed for Chronic Wasting Disease (CWD) infections in deer and elk facing predation by cougars and grey wolves in the Greater Yellowstone Ecosystem. The results indicate that predators can reduce CWD outbreak size, especially if selecting for infected individuals. CWD is an always fatal disease and this work suggests that predators could reduce disease burdens in cervids. The model is also applicable to other systems and promises to further our understanding of the role of predation on disease in prey, as well as drive future empirical studies.

Genetic Characterization of Microsporum canis Clinical Isolates in the United States.

Microsporum canis is the primary agent causing dermatophytosis in cats, and also infects humans, dogs, and other species. Assessment of genetic variation among M. canis isolates in the United States has not been conducted. Further, M. canis mating type and assessment of disease severity associated with genotypic characteristics have not been rigorously evaluated. We therefore isolated M. canis from 191 domestic cats across the US and characterized genotypes by evaluation of ITS sequence, MAT locus, and microsatellite loci analysis. The genes SSU1 and SUB3, which are associated with keratin adhesion and digestion, were sequenced from a subset of isolates to evaluate potential genetic associations with virulence. Analysis of microsatellite makers revealed three M. canis genetic clusters. Both clinic location and disease severity were significant predictors of microsatellite variants. 100% of the M. canis isolates were MAT1-1 mating gene type, indicating that MAT1-2 is very rare or extinct in the US and that asexual reproduction is the dominant form of replication. No genetic variation at SSU1 and SUB3 was observed. These findings pave the way for novel testing modalities for M. canis and provide insights about transmission and ecology of this ubiquitous and relatively uncharacterized agent.

Multi-population puma connectivity could restore genomic diversity to at-risk coastal populations in California.

Urbanization is decreasing wildlife habitat and connectivity worldwide, including for apex predators, such as the puma (Puma concolor). Puma populations along California's central and southern coastal habitats have experienced rapid fragmentation from development, leading to calls for demographic and genetic management. To address urgent conservation genomic concerns, we used double-digest restriction-site associated DNA (ddRAD) sequencing to analyze 16,285 genome-wide single-nucleotide polymorphisms (SNPs) from 401 pumas sampled broadly across the state. Our analyses indicated support for 4-10 geographically nested, broad- to fine-scale genetic clusters. At the broadest scale, the four genetic clusters had high genetic diversity and exhibited low linkage disequilibrium, indicating that pumas have retained genomic diversity statewide. However, multiple lines of evidence indicated substructure, including 10 finer-scale genetic clusters, some of which exhibited fixed alleles and linkage disequilibrium. Fragmented populations along the Southern Coast and Central Coast had particularly low genetic diversity and strong linkage disequilibrium, indicating genetic drift and close inbreeding. Our results demonstrate that genetically at risk populations are typically nested within a broader-scale group of interconnected populations that collectively retain high genetic diversity and heterogenous fixations. Thus, extant variation at the broader scale has potential to restore diversity to local populations if management actions can enhance vital gene flow and recombine locally sequestered genetic diversity. These state- and genome-wide results are critically important for science-based conservation and management practices. Our nested population genomic analysis highlights the information that can be gained from population genomic studies aiming to provide guidance for the conservation of fragmented populations.

Hunting alters viral transmission and evolution in a large carnivore.

Hunting can fundamentally alter wildlife population dynamics but the consequences of hunting on pathogen transmission and evolution remain poorly understood. Here, we present a study that leverages a unique landscape-scale quasi-experiment coupled with pathogen-transmission tracing, network simulation and phylodynamics to provide insights into how hunting shapes feline immunodeficiency virus (FIV) dynamics in puma (Puma concolor). We show that removing hunting pressure enhances the role of males in transmission, increases the viral population growth rate and increases the role of evolutionary forces on the pathogen compared to when hunting was reinstated. Changes in transmission observed with the removal of hunting could be linked to short-term social changes while the male puma population increased. These findings are supported through comparison with a region with stable hunting management over the same time period. This study shows that routine wildlife management can have impacts on pathogen transmission and evolution not previously considered.

Parasites as conservation tools.

Parasite success typically depends on a close relationship with one or more hosts; therefore, attributes of parasitic infection have the potential to provide indirect details of host natural history and are biologically relevant to animal conservation. Characterization of parasite infections has been useful in delineating host populations and has served as a proxy for assessment of environmental quality. In other cases, the utility of parasites is just being explored, for example, as indicators of host connectivity. Innovative studies of parasite biology can provide information to manage major conservation threats by using parasite assemblage, prevalence, or genetic data to provide insights into the host. Overexploitation, habitat loss and fragmentation, invasive species, and climate change are major threats to animal conservation, and all of these can be informed by parasites.

Los Parásitos como Herramienta de Conservación Resumen El éxito de los parásitos depende típicamente de la relación cercana con uno o más hospederos; por lo tanto, las características de la infección parasitaria tienen potencial para proporcionar detalles indirectos de la historia natural del hospedero y son biológicamente relevantes para la conservación animal. La caracterización de las infecciones parasitarias ha sido útil para definir a las poblaciones hospederas y ha servido como sustituto para la evaluación de la calidad ambiental. Los estudios innovadores de la biología de parásitos pueden proporcionar información para manejar las principales amenazas a la conservación mediante la información proporcionada por el conjunto de parásitos, su prevalencia o genética que proporciona conocimiento sobre el hospedero. La sobreexplotación, la pérdida del hábitat y la fragmentación, las especies invasoras y el cambio climático son las principales amenazas para la conservación animal y a todas pueden ser informadas mediante los parásitos.

INVESTIGATING ASSOCIATIONS AMONG RELATEDNESS, GENETIC DIVERSITY, AND CAUSES OF MORTALITY IN SOUTHERN SEA OTTERS (ENHYDRA LUTRIS NEREIS).

Southern sea otter (Enhydra lutris nereis) population recovery is influenced by a variety of factors, including predation, biotoxin exposure, infectious disease, oil spills, habitat degradation, and resource limitation. This population has also experienced a significant genetic bottleneck, resulting in low genetic diversity. We investigated how two metrics, familial relatedness and genetic diversity, are correlated with common causes of mortality in southern sea otters, including cardiomyopathy, acanthocephalan (Profilicollis spp.) peritonitis, systemic protozoal infection (Toxoplasma gondii and Sarcocystis neurona), domoic acid intoxication, end-lactation syndrome, and shark bite. Microsatellite genetic markers were used to examine this association in 356 southern sea otters necropsied from 1998 to 2012. Significant associations with genetic diversity or familial relatedness (P<0.05) were observed for cardiomyopathy, acanthocephalan peritonitis, and sarcocystosis, and these associations varied by sex. Adult male cardiomyopathy cases (n=86) were more related than the null expectation (P<0.049). Conversely, female acanthocephalan peritonitis controls (n=110) were more related than the null expectation (P<0.004). Including genetic diversity as a predictor for fatal acanthocephalan peritonitis in the multivariate logistic model significantly improved model fit; lower genetic diversity was associated with reduced odds of sea otter death due to acanthocephalan peritonitis. Finally, male sarcocystosis controls (n=158) were more related than the null expectation (P<0.011). Including genetic diversity in the multivariate logistic model for fatal S. neurona infection improved model fit; lower genetic diversity was associated with increased odds of sea otter death due to S. neurona. Our study suggests that genetic diversity and familial relatedness, in conjunction with other factors such as age and sex, may influence outcome (survival or death) in relation to several common southern sea otter diseases. Our findings can inform policy for conservation management, such as potential reintroduction efforts, as part of species recovery.

SARS-CoV-2 evolution in animals suggests mechanisms for rapid variant selection.

SARS-CoV-2 spillback from humans into domestic and wild animals has been well documented, and an accumulating number of studies illustrate that human-to-animal transmission is widespread in cats, mink, deer, and other species. Experimental inoculations of cats, mink, and ferrets have perpetuated transmission cycles. We sequenced full genomes of Vero cell-expanded SARS-CoV-2 inoculum and viruses recovered from cats (n = 6), dogs (n = 3), hamsters (n = 3), and a ferret (n = 1) following experimental exposure. Five nonsynonymous changes relative to the USA-WA1/2020 prototype strain were near fixation in the stock used for inoculation but had reverted to wild-type sequences at these sites in dogs, cats, and hamsters within 1- to 3-d postexposure. A total of 14 emergent variants (six in nonstructural genes, six in spike, and one each in orf8 and nucleocapsid) were detected in viruses recovered from animals. This included substitutions in spike residues H69, N501, and D614, which also vary in human lineages of concern. Even though a live virus was not cultured from dogs, substitutions in replicase genes were detected in amplified sequences. The rapid selection of SARS-CoV-2 variants in vitro and in vivo reveals residues with functional significance during host switching. These observations also illustrate the potential for spillback from animal hosts to accelerate the evolution of new viral lineages, findings of particular concern for dogs and cats living in households with COVID-19 patients. More generally, this glimpse into viral host switching reveals the unrealized rapidity and plasticity of viral evolution in experimental animal model systems.

MrIML: Multi-response interpretable machine learning to model genomic landscapes.

We introduce a new R package "MrIML" ("Mister iml"; Multi-response Interpretable Machine Learning). MrIML provides a powerful and interpretable framework that enables users to harness recent advances in machine learning to quantify multilocus genomic relationships, to identify loci of interest for future landscape genetics studies, and to gain new insights into adaptation across environmental gradients. Relationships between genetic variation and environment are often nonlinear and interactive; these characteristics have been challenging to address using traditional landscape genetic approaches. Our package helps capture this complexity and offers functions that fit and interpret a wide range of highly flexible models that are routinely used for single-locus landscape genetics studies but are rarely extended to estimate response functions for multiple loci. To demonstrate the package's broad functionality, we test its ability to recover landscape relationships from simulated genomic data. We also apply the package to two empirical case studies. In the first, we model genetic variation of North American balsam poplar (Populus balsamifera, Salicaceae) populations across environmental gradients. In the second case study, we recover the landscape and host drivers of feline immunodeficiency virus genetic variation in bobcats (Lynx rufus). The ability to model thousands of loci collectively and compare models from linear regression to extreme gradient boosting, within the same analytical framework, has the potential to be transformative. The MrIML framework is also extendable and not limited to modelling genetic variation; for example, it can quantify the environmental drivers of microbiomes and coinfection dynamics.

Feline Leukemia Virus (FeLV) Endogenous and Exogenous Recombination Events Result in Multiple FeLV-B Subtypes during Natural Infection.

Feline leukemia virus (FeLV) is associated with a range of clinical signs in felid species. Differences in disease processes are closely related to genetic variation in the envelope (env) region of the genome of six defined subgroups. The primary hosts of FeLV are domestic cats of the Felis genus that also harbor endogenous FeLV (enFeLV) elements stably integrated in their genomes. EnFeLV elements display 86% nucleotide identity to exogenous, horizontally transmitted FeLV (FeLV-A). Variation between enFeLV and FeLV-A is primarily in the long terminal repeat (LTR) and env regions, which potentiates generation of the FeLV-B recombinant subgroup during natural infection. The aim of this study was to examine recombination behavior of exogenous FeLV (exFeLV) and enFeLV in a natural FeLV epizootic. We previously described that of 65 individuals in a closed colony, 32 had productive FeLV-A infection, and 22 of these individuals had detectable circulating FeLV-B. We cloned and sequenced the env gene of FeLV-B, FeLV-A, and enFeLV spanning known recombination breakpoints and examined between 1 and 13 clones in 22 animals with FeLV-B to assess sequence diversity and recombination breakpoints. Our analysis revealed that FeLV-A sequences circulating in the population, as well as enFeLV env sequences, are highly conserved. We documented many recombination breakpoints resulting in the production of unique FeLV-B genotypes. More than half of the cats harbored more than one FeLV-B variant, suggesting multiple recombination events between enFeLV and FeLV-A. We concluded that FeLV-B was predominantly generated de novo within each host, although we could not definitively rule out horizontal transmission, as nearly all cats harbored FeLV-B sequences that were genetically highly similar to those identified in other individuals. This work represents a comprehensive analysis of endogenous-exogenous retroviral interactions with important insights into host-virus interactions that underlie disease pathogenesis in a natural setting. IMPORTANCE Feline leukemia virus (FeLV) is a felid retrovirus with a variety of disease outcomes. Exogenous FeLV-A is the virus subgroup almost exclusively transmitted between cats. Recombination between FeLV-A and endogenous FeLV analogues in the cat genome may result in emergence of largely replication-defective but highly virulent subgroups. FeLV-B is formed when the 3' envelope (env) region of endogenous FeLV (enFeLV) recombines with that of the exogenous FeLV (exFeLV) during viral reverse transcription and integration. Both domestic cats and wild relatives of the Felis genus harbor enFeLV, which has been shown to limit FeLV-A disease outcome. However, enFeLV also contributes genetic material to the recombinant FeLV-B subgroup. This study evaluates endogenous-exogenous recombination outcomes in a naturally infected closed colony of cats to determine mechanisms and risk of endogenous retroviral recombination during exogenous virus exposure that leads to enhanced virulence. While FeLV-A and enFeLV env regions were highly conserved from cat to cat, nearly all individuals with emergent FeLV-B had unique combinations of genotypes, representative of a wide range of recombination sites within env. The findings provide insight into unique recombination patterns for emergence of new pathogens and can be related to similar viruses across species.

SARS-CoV-2 evolution in animals suggests mechanisms for rapid variant selection.

SARS-CoV-2 spillback from humans into domestic and wild animals has been well-documented. We compared variants of cell culture-expanded SARS-CoV-2 inoculum and virus recovered from four species following experimental exposure. Five nonsynonymous changes in nsp12, S, N and M genes were near fixation in the inoculum, but reverted to wild-type sequences in RNA recovered from dogs, cats and hamsters within 1-3 days post-exposure. Fourteen emergent variants were detected in viruses recovered from animals, including substitutions at spike positions H69, N501, and D614, which also vary in human lineages of concern. The rapidity of in vitro and in vivo SARS-CoV-2 selection reveals residues with functional significance during host-switching, illustrating the potential for spillback reservoir hosts to accelerate evolution, and demonstrating plasticity of viral adaptation in animal models. ONE-SENTENCE SUMMARY: SARS-CoV-2 variants rapidly arise in non-human hosts, revealing viral evolution and potential risk for human reinfection.

A natural laboratory to elucidate the evolution of endogenous-exogenous retroviral interactions.

The advent of whole genome sequencing has revealed much about the genomes of animals including the relatively large percentage of the genome consisting of endogenous retroviruses (ERV; International Human Genome Sequencing Consortium, 2001). An ERV arises when a retrovirus integrates into a host germ cell genome through normal infection processes. Germline infections can be transmitted to offspring through Mendelian inheritance and at times become fixed elements of the host genome (Weiss, 2006). At their inception, these endogenized retroviruses maintain all the functions of their exogenous progenitors and can produce infectious virus (Figure 1). Mutations in the ERV randomly accumulate over time and can lead to the loss of the deleterious effects. In addition, the ERV can provide benefits to the host often via limiting exogenous viral infections (Chiu and VandeWoude, 2020a). However, most of these endogenous viruses are evolutionary relics representing historical infections and have no contemporary exogenous virus. This limits the opportunities to understand the evolution of ERVs and their interactions with exogenous viruses. In this issue of Molecular Ecology, Quigley et al. (2020) use a novel approach to show that koala retrovirus (KoRV) is undergoing endogenization along a geographic gradient with a variety of exogenous variants dispersed across the landscape. This system provides an opportunity to further elucidate the complex mechanisms in which endogenous and exogenous viruses interact and follow the evolution of an ERV in real time.

Host relatedness and landscape connectivity shape pathogen spread in the puma, a large secretive carnivore.

Urban expansion can fundamentally alter wildlife movement and gene flow, but how urbanization alters pathogen spread is poorly understood. Here, we combine high resolution host and viral genomic data with landscape variables to examine the context of viral spread in puma (Puma concolor) from two contrasting regions: one bounded by the wildland urban interface (WUI) and one unbounded with minimal anthropogenic development (UB). We found landscape variables and host gene flow explained significant amounts of variation of feline immunodeficiency virus (FIV) spread in the WUI, but not in the unbounded region. The most important predictors of viral spread also differed; host spatial proximity, host relatedness, and mountain ranges played a role in FIV spread in the WUI, whereas roads might have facilitated viral spread in the unbounded region. Our research demonstrates how anthropogenic landscapes can alter pathogen spread, providing a more nuanced understanding of host-pathogen relationships to inform disease ecology in free-ranging species.